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Abstract. We present Far Ultraviolet Spectroscopic Explorer observations of the star-forming galaxy Mrk 54 at 
z = 0.0448. The Lyman continuum radiation is not detected above the Hi absorption edge in our Galaxy. An 
upper limit is evaluated by comparison with the background measured in regions of the detector adjacent to the 
observed spectrum. A spectral window of 16 A, reasonably free of additional Hi Lyman series line absorption 
is used. No correction is needed for molecular hydrogen absorption in our Galaxy but a foreground extinction 



of 0.29 mag is accounted for. An upper limit of 6.15 x 10 erg cm 



A is obtained for the flux at 



900 A in the rest frame of Mrk 54. By comparison with the number of ionizing photons derived from the Ha 
flux, this limit translates into an upper limit of / csc < 0.062 for the fraction of Lyman continuum photons that 
escape the galaxy without being absorbed by interstellar material. This limit compares with the limits obtained 
in three other nearby galaxies and is compatible with the escape fractions predicted by models. The upper limits 
obtained in nearby galaxies contrasts with the detection of Lyman continuum flux in the composite spectrum of 
Lyman-break galaxies at z w 3.4. The difficulties and implications of a comparison are discussed. 
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1. Introduction 

It is not yet clear whether hot and massive stars forming 
in galaxies contribute significantly to the ionizing back- 
ground radiation in the universe, and how this contribu- 
tion evolves as a function of reds hift ( e.g. Bechtold et al. 
1981 Miralda-Escude fc Os triker |l99dj , Meiksi n fc M adau 
19931 , Mada u & Shull |l996[ Haardt & Madau |1996[ Shull 
et al. 1999). Specifically, at redshifts z > 3, the early for- 



( 1997| ) obtained upper limits on the Lyman continuum 
(LyC) radiation in four nearby star-forming galaxies. By 
comparison with the Ha emission these data were inter- 
preted in terms of limits on the LyC escape fraction, a 
parameter giving the fraction of hydrogen-ionizing pho- 
tons effectively released into the IGM. Recently Steidel 
et al. (2001) reported the detection of the LyC radia- 



mation of galaxies is expected to compensate for the de- 
cline in the quasar contribution and to play a role in the 
re-ionization of the IGM (e.g. Madau et al. 1999] ). 

Direct observations of galaxies below the Lyman break 
have been scarce so far. With the Hopkins Ultraviolet 
Telescope (HUT), Leitherer et al. ( |1995| ), Hurwitz et al. 



tion in a composite spectrum of Lyman break galaxies 
at z ~ 3.4 that was also interpreted in terms of LyC es- 
cape fraction but by comparison with the (1500 A) UV 
continuum. A preliminary account of the Space Telescope 
Imaging Spectrograph (STIS) observations of galaxies in 
the Hubble Deep Field (HDF), with implications on the 
LyC escape fraction, has been presented by Ferguson 
(2001) while this paper was nearing completion. 
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In the absence of a large number of observations that 
would directly provide the LyC luminosity function and 
the LyC luminosity density of galaxies, the LyC escape 
fraction is seen as a crucial parameter. Combined with 
the Ha or UV luminosity densities of galaxies or with 
the stellar ionizing radiation calculated from evolutionary 
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synthesis models (e.g. Bruzual & Chariot 1993, Leitherer 



et al. 1999) or from the rate of chemical enrichment in 
the universe (Cowie 1988, Songaila et al. 199C, Madau & 
Shull 1996Q , the LyC escape fraction provides the amount 
of ionizing radiation effectively released into the IGM by 
the galaxies. As this parameter encapsulates a number of 
complex and random factors it is probably highly variable 
from galaxy to galaxy; its full understanding and sound 
utilisation would also require a large number of observa- 
tions. 

The Far Ultraviolet Spectroscopic Explorer (FUSE) has 
recently opened again an access into the far ultraviolet 
down to 905 A, allowing the possibility to observe the LyC 
radiation of low redshift galaxies above the 912 A Lyman 
limit of Hi photoelectric absorption in our Galaxy. As 
sho wn by the re-analysis of HUT data by Hurwitz et 
al. ( 1997 ), getting rid of photoelectric absorption in our 
Galaxy above the Lyman edge is not enough and gas-phase 
absorption in our Galaxy, essentially from convergent H I 
Lyman series and the Lyman and Werner bands of molec- 
ular hydrogen has to be accounted for in addition to dust 
extinction. It is an area where the high spectral resolu- 
tion of FUSE can help even though the increased back- 
ground contribution per wavelength unit makes the de- 
tection of faint continua more difficult than at low resolu- 
tion. In the following we report the FUSE observations of 
a star-forming galaxy Mrk 54, at a redshift z = 0.0448 that 
places the Lyman break at 952.5 A above the Galactic 
LyS absorption feature, reasonably beyond the convergent 
H I Lyman series. 

2. Observations and Data Processing 

The observations of Mrk 54 (GI program A052) were ob- 
tained on Feb 19 - 20, 2000 in the time-tagged photon 
address mode with the object in the large (30" x 30") 
aperture. The observations were split into 15 exposures be- 
tween occultation periods and passages through the South 
Atlantic Anomaly, and generally extended over spacecraft 
night and day. The total duration was 27502 seconds. 
Details of the FUSE instrument and on-orbit performance 
have b een given by Moos et al. ( [2000 ) and Sahnow et al. 
( 120001) . 

The summed and calibrated spectrum of Mrk 54 is 
consistent at long wavelength with the IUE observations 
(Kinney et al. 1993| ) and reaches a maximum flux of about 

3. x 10~ 14 erg cm" 2 s" 1 A" 1 at ~ 1100 A. At wavelengths 
shorter than 1000 A the flux decreases rapidly and be- 
low 960 A the signal, if any, becomes extremely weak 
(Figure 1); actually it drops to below zero between the 
Lyman-series airglow lines with the uniform background 
correction of 1 count cm -2 s _1 applied by the standard 
calibration pipeline. The faint level between the redshifted 
Lyman break of the galaxy 911.7 (1+0.0448) = 952.5 
A and about 960 A is probably due to the accumula- 
tion of Lyman series absorption lines in the galaxy itself. 
These lines are resolved from Ly/3 to Ly7 (Ly7 and Ly/3 
are not shown in Fig 1). 
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Fig. 1. FUSE spectrum of Mrk 54 resulting from the co- 
addition of 15 SiC2A spectra for a total exposure time of 
27502 s and at 0.2 A linear rebinning. The Lyman series 
absorption lines resolved in Mrk 54 are marked except Ly6 
confused with the Ly7 airglow line at 972.5 A. 
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Fig. 2. Details of the previous spectrum at 0.1 A rebin- 
ning, showing (dotted line) the three spectral windows (16 
A total) used for obtaining an upper limit to the LyC of 
Mrk 54. Tentative i dentifi cations of the airglow lines based 
on Feldman et al. (p000j) are indicated. 



For the purpose of discussing the LyC radiation of the 
galaxy Mrk 54 we have concentrated our analysis on the 
SiC channels and specifically the SiClB and SiC2A spec- 
tra. We have screened the raw ttag data of segments IB 
and 2A for the presence of so-called bursts, using for each 
of the 15 exposures both the images of the events and 
the time series of the total number of counts. The latter 
has shown some increase of the count rate at the end of 
a few exposures: as these features were not typical burst 
events we decided not to remove them from the ttag data 
(even for the most significant in exposure 13). We made 
no attempt at isolating night-only data. 

In the absence of any clear signal from the redshifted 
LyC of the galaxy (shortward of 952.5 A) we have tried 
to place an upper limit to this flux. A first step is to de- 
termine spectral windows that are free from additional 
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absorption by neutral hydrogen in our Galaxy (Hurwitz 
et al. 1997). As discussed by these latter authors, it is a 
complex issue given the lack of information regarding low 
column density clouds at velocities departing significantly 
from that of the bulk of the Hi. The examples given by 
Hurwitz et al. (1997) and by Lockman & Savage (1995) as 
well as the examination of the list of high-velocity clouds 



of Stark et al. (1992) lead us to conclude that removing ± 
200 km s' 1 on either side of the Lyman series absorption 
lines is probably safe enough. This selection removes the 
H I nightglow lines sitting right at the rest wavelengths of 
the potential absorption lines. A comparison of our spec- 
tru m wit h the reference airglow spectrum of Feldman et 
al. ( EOQCI) shows that we should also avoid a few O I lines 
close to H I lines and discard a slightly larger wavelength 
domain shortward of each Lyman series absorption lines. 



This corresponds to a velocity of 



500 km s , mak- 



ing our previous limit even safer on the side of negative 
velocities. With the constraints adopted it is not possi- 
ble to extend our study shortward of Ly7 (926.2 A). We 
practically end up with a window of 16 A total, split into 
three domains, 927 - 929 A, 931.5 - 936 A and 938.5 - 
948 A by the Lyman series lines Ly6 930.7 A, Lye 937.8 
A and Ly<5 949.7 A (Figure 2). We have calculated the 
total raw counts in this window for each of the 15 SiClB 
and 15 SiC2A extracted spectra. 

These raw counts can be compared with those obtained 
in reference zones, defined in the 15 images of segments 
IB and 2A built from the corresponding raw ttag data 
(Figure 3). These reference zones, supposed to represent 
the background at the time of observations, have been 
searched as close as possible from the detector areas where 
the sky observations are made. Our search was guided by 
the examination of profiles integrated over a large num- 
ber of lines or columns of the co-added (15) images in 
order to avoid detector edge effects and airglow lines in 
the smaller apertures. We end up with three areas of re- 
spectively 542101, 1023701 and 737751 pixels in the IB 
segment and 322361, 610151 and 443631 pixels in the 2A 
segment. In these zones the counts are found to follow 
Poisson statistics extremely well. 

Three features stand out from Figure 3. First, given 
their (1 a) error bars the three background measurements 
are neither significatively different in each exposures nor 
show systematic variation pattern. This is taken as a good 
indication that the background is reasonably uniform at 
least at the scale and in the zones we are using. In contrast 
there are significant variations of the background measure- 
ments (departing from purely Poisson statistics) from ex- 
posure to exposure that can be understood by difference 
of observing conditions. Third, the target measurements 
follow the same trend as the background from exposure to 
exposure (except may be for exposure 7) and, given their 
larger error bars, do not appear significantly different from 
the background measurements. 

These findings justify the use of our evaluation of the 
background in each image for correcting the sky measure- 
ment. The error bars on the resulting net counts in each 
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Fig. 3. Comparison of the sky raw counts (open circle) in 
the 16 A window defined in the text with the background 
counts in three reference zones (solid dots). These zones 
are displayed in the same relative order for each 15 expo- 
sures. For the purpose of comparison all counts have been 
expressed into units of counts cm' 2 s" 1 using the respec- 
tive exposure times, the surface of the extracting windows 
for the sky counts, the number of pixels and their linear 
size (respectively 5.98 x 9.08 /im for the IB segment and 
5.96 x 14.75 /zm for the 2A segment) for the background 
zones, la Poisson error bars are given. 



exposure combine quadratically the error bar on the sky 
measurements and those on the background scaled to the 
surface of the sky measurements; they are dominated by 
the uncertainties on the sky measurements based on a 
smaller number of pixels. The error bars of each 15 ex- 
posures are then combined quadratically and divided by 
the total exposure time to give the dispersion on the net 
count rate. We obtain a mean net count rate of 1.32 x 
10 -3 count s _1 with a 1 a dispersion of 1.40 x 10~ 3 count 
s _1 and —1.13 x 10 -3 count s _1 with a 1 a dispersion of 
1.43 x 10~ 3 count s -1 in the 16 A window of segments IB 
and 2A respectively. The fact that the mean count rates 
are smaller than the dispersion confirms quantitatively the 
trends discussed above. In these conditions the net count 
rate dispersion can be used to set an upper limit on the 
flux of Mrk 54 at w 900A (rest-frame). Combining the 
two segments and converting the counts into flux (using 
an effective area of the order of 17 cm 2 at ~ 940 A ac- 
cording to on-orbit performance reported by Sahnow et al. 



[2U00 ) we get a 3 a upper limit of 4.7 x 10 16 erg cm 2 s 1 



A -1 . We have carried out all these calculations in counts 
rather than calibrated flux units in order to stay closer 
to basic count statistics. This advantage offsets the slight 
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inaccuracy of converting counts into fluxes for the entire 
spectral window. 



3. Analysis and discussion 

3.1. Correction to the upper limit 

We need first to correct our upper limit flux from any 
absorption that would not take place in the object but 
along its line of sight in the intergalactic medium or in 
our Galaxy. At the low redshift of Mrk 54 the intergalac- 
tic medium opacity can be neglected. As to the gas-phase 
absorption in our Galaxy, we have seen that the choice of 
our 16 A spectral window makes significant absorption by 
neutral hydrogen unlikely. Most of the known interstellar 
metal lines, especially O I, are also found to be avoided and 
the equivalent width of the corresponding absorption is 
negligible with respect to the width of the window. For the 
molecular hydrogen that FUSE has shown to be present 
along most of the extragalactic lines of sight (Shull et al. 
2000| ) we have searched for the presence of the most signif- 
icant R(0) and R(l) lines of the low-rotational levels of the 
Lyman and Werner bands, following the ide ntific ations in 
the line of sight of ESO 14 1-055 (Shull et al. fc000| ) and the 
tables of Barnstedt et al. ( |2000|) . We found no clear iden- 
tifications even in the brightest parts of the spectrum. By 
comparison with the e xamples of ESO 141-055 and Mrk 
876 (Shull et al. |2000|) we conclude that N(R 2 ) < 1 x 



10 cm in the direction of Mrk 54. A difficulty of this 
comparison comes from the fact that the spectra of ESO 
141-055 and Mrk 876 have better signal to noise ratio than 
Mrk 54 (the objects are approximately two times brighter 
and the exposure times two times longer). We have used 
the on-line spectral simulator to explore how the detec- 
tion is affected by these differences and found that a H2 
column density of 1 x 10 18 cm~ 2 would still have been 
detected in the conditions of our Mrk 54 spectrum. If we 
give up any comparison with observed spectra and use the 
simulator alone a lower limit of 1 x 10 17 cm -2 is found. 
We have then used the simulator to directly calculate the 
fraction of energy that such H2 column densities would ab- 
sorb into our specific window of 16 A; we found a fraction 
of 0.015 for 1 x 10 17 cm~ 2 and 0.043 for the more conser- 
vative limit of 1 x 10 18 cm -2 . We have therefore applied 
no correction for H2 absorption to our upper limit on the 
LyC of Mrk 54. This case may well be a fortunate circum- 
stance given the galactic H2 column densities cur rently 
reported for extraga lactic lines of sight (Shull ct al. 200C , 
Vidal-Madjar et al. 2000 ). In the absence of a control of 
the amount of H2 along the line of sight, that FUSE now 
offers, previous measurements may well have been affected 
by this problem. If accounted for, the H2 contamination 
would have led to even less restrictive limits on th e LyC 
escape fraction than established by Hurwitz et al. ( 1997 ). 

Last, the Galactic dust extinction is a more severe and 
more uncertain factor. In contrast to the current thinking 
for the past 17 years, it is very likely that a residual ex- 
tinction be present at high Galactic latitudes, even in the 



directions with the lowest H 1 column densities (Schlegel et 
al. [1998D - For Mrk 54 we have adopted E{B - V) = 0.015 
as given by the NASA Extragalactic Database on the basis 
of the latter reference. This color excess is consistent with 
a low fraction of hydrogen in the molecular state (Savage 
et al. 1977) and therefore with previous limits placed on 



the H2 column density along the line of sight to Mrk 54. 
This results in a foreground extinction of 0.29 mag in our 
far-UV window (centered at ~ 940 A) using an extrapo- 
lation of the parameterized extinction law of Cardelli et 
al. (1989) shortward of 1000 A. Such an extra polati on is 
supported by the measurements of Buss ct al. ( 1994 ). We 
have therefore to account for a factor 1.306 of absorption 
that does not occur in Mrk 54; our upper limit on the 
/(900) flux of Mrk 54 is increased to 6.15 x lO" 16 erg 
cm~ 2 s _1 A -1 . 



3.2. Evaluation of the LyC escape fraction 

We can now proceed with the calculation of an upper limit 
on the LyC escape fraction / esc defined as the fraction of 
emitted 900 A photons that escapes the galaxy without 
being absorbed by interstellar material. We need first to 
have a relation between the /(900) flux and the total num- 
ber of LyC photons since it is this latter quantity that 
can be physical ly con strained by the observed Ha flux. 
Leitherer et al. ( 1995 ) have shown that there is a narrow 
relation relatively independent of star formation histories 
and initial mass functions between the luminosity L(900) 
at about 900 A and the total number A^y of LyC pho- 
tons of a burst population: log(7V Ly /i(900)) = 13.28±0.16 
(photons A erg" 1 ). 

Assuming 0.45 as the number of Ha photons per re- 
combination (case B recombination at 10 4 K), we write the 
number of LyC photons A^y, corrected by the fraction of 
those escaping photoelectric absorption (/l) or trapped 
by dust before ionization (f2), 

fl - /2) = 7.34 x 10 11 47rD 2 /(Ha)10 a4 ' 4(HQ) 



^Ly(l 

with /(Ha) the observed Ha flux (in erg cm -2 s _1 ), 
A(Ha) the extinction at Ha in the observed galaxy and 
D the distance (in cm). The luminosity L(900) is 

i(900)/ csc = 47r£ 2 /(900) 

The continuum flux /(900) is corrected for gas-phase ab- 
sorption in the Galaxy and foreground dust extinction. It 
should be noted that / osc may be different from fl because 
of photons escaping photoelectric absorption but absorbed 
by dust and the frequency dependence of photoelectric ab- 
sorption. Combination of above equations gives: 

26/(900)10-° 4 ' 4 ( Hq )(1 - /I - /2) 



fc: 



(1) 



/(Ha) 

As we are dealing with an upper limit calculation in the 
case of Mrk 54 and we expect / eS c, /l, and /2 to be small 
it makes sense to write 

26/(900) l0~oAA(H a ) 



/esc 



/(Ha) 



(2) 
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An Ha flux of 2.6 x 10~ 13 erg cm~ 2 s -1 has been ob- 
tained for Mrk 54 from spectrophotometric observations 
(Bosclli 200C). An extinction A(Ha) = (although a value 
of 0.35 was obtained from the Balmer decrement measured 
along one slit position) has been adopted in order to ac- 
count for the uncertainty on the Ha flux in the upper limit 
calculation of equation (2). We finally obtain / esc < 0.062 
for Mrk 54. 

This result adds another significant limit to three other 
cases with limits of 3.2%, 5.2% and 11% as obtained with 



the interstellar media of starburst galaxies (e.g. Heckman 
000 and references therein) and play a role in the escape 
of Lya emission (e.g. Kunth et al. 1998). As the Lya emis- 



HUT and the re-analysis by Hurwitz et al. ( |1997| ) (Mrk 66 
with a limit of only 57% is left over from this comparison). 
All the values obtained so far are fully compatible with 
current estimates in the range 2% - 10% based either on 



theoretical models (Dove & Shull 1994, Dove et al. 2000) 



or the implications of Ha observatio ns of the Magellanic 
Stream (Bland-H awtho rn & Maloney |l999|) or NGC 3067 
(Tumlinson et al. |l999| ). 



3.3. Neutral hydrogen 

Neutral hydrogen is expected in galaxies where active star 
formation is taking place and its impact on the escape of 
ionizing photons is thought to depend heavily on topol- 
ogy. In Mrk 54 the neutral hydrogen is directly seen along 
the line of sight as Lyman series absorption lines (down 
to Ly7) at the redshift of the parent galaxy. These fea- 
tures are accompanied by a number of metal absorption- 
lines with different possible origins (interstellar gas, stel- 
lar photospheres and stellar winds) as currently found in 
starburst galaxies (e.g. Heckman et al. 199S and references 
therein). Most of the analyses so far have concentrated on 
the features (including Lya) accessible with IUE or the 
spectrographs of the HST. To the best of our knowledge, 
the discussion of features blueward of Lya have been lim- 
ited to the O VI + Ly/3 + C n profile in the four starbursts 
observed with HUT (Gonzalez-Delgado et al. |1997[ |199Sp . 

Leaving a full discussion of absorption features in Mrk 
54 to another context, we have only tried to interpret the 
Hi Lyman series absorption profiles. Within the uncer- 
tainties of background subtraction, the lines appear al- 
most black at their centers (this is less clear with the Ly(3 
and Ly7 lines than with the three lines displayed in Figure 
1) and indicate a large covering fraction of neutral hydro- 
gen. The determination of the neutral hydrogen column 
density implied by the absorption profiles is however diffi- 
cult because of the low signal-to-noise ratio and the possi- 
ble blending with galactic interstellar lines (not to speak of 
nightglow emission lines) resulting in very uncertain con- 
tinuum levels: the profile fitting allows column densities 
of the order of 10 21 cm -2 with low b Doppler parameter 
as well as low column densities with b up to ~ 300 km 
s _1 . Large velocity spreads have already been reported in 
nearby and high redshift starburst galaxies (e.g. Gonzalez- 
Delgado et al. |199S| , Pettini et al. |2000| ). 

This finally leads us to place the discussion in the con- 
text of the large-scale outflows that have been revealed in 



sion often appears redshifted with respect to absorption 
features the simplest interpretation is that the only Lya 
photons that escape unabsorbed are those backscattered 
from the far side of the outflow whereas the approaching 
part is seen in absorption against the stellar continuum 
(Pettini et al. 200C). In this picture the LyC photons can- 



not escape in the same way as the backscattered Lya pho- 
tons but it is likely that the galactic winds generate holes 
in the H I distribution through which LyC photons can 
also escape in different directions. A more complete model, 
including time evolution effect and accounting for the va- 
riety of Lya absorp tion and emission profiles observed by 
Kunth et al. ( |l998 ), has been developed by Tenorio-Tagle 
et al. ( 1999] ); in the phase where the conical Hn region 
extends to the galaxy outer edge, LyC and Lya photons 
should escape in the same direction. With only upper lim- 
its as obtained so far, and only five nearby galaxies, it is 
difficult to constrain one particular model by examining 
how the LyC escape fraction correlates with the Lya es- 
cape and the blend of absorption and emission at Lya; for 
instance the two objects with the tightest LyC escape frac- 
tion upper limit (Mrk 496 and IRAS 08339) have both Lya 
in emission. The only sure conclusion of this discussion is 
the role of anisotropics that increase the randomness of 
the LyC escape fraction for the observer. 

4. Comparison with the Lyman break galaxies 

The upper limits (< 10%) obtained at z rj are in stark 
contrast with a LyC escape fraction larger than 50% re- 
ported by Steidel et al. (2001) from the detection of the 
LyC radiation in a composite spectrum of Lyman break 
galaxies at z — 3.4. In Steidel et al. the LyC escape fraction 
is normalized by the fraction of 1500 A photons that es- 
capes. The surprisingly large value (> 50%) results from 
the fact that the /(1500)//(900) flux ratio is, after the 
necessary correction for the intergalactic medium opacity 
at high redshifts, close to that predicted from spectral syn- 
thesis models without any LyC self-absorption from neu- 
tral hydrogen in the galaxy. The escape fractions used at 
low and high redshifts are therefore different both in ref- 
erence wavelength (the escape fraction for nearby galaxies 
refers to Ha vs. 1500 A for the Lyman break galaxies) 
and in the treatment of dust extinction (dust extinction 
is accounted for at Ha in nearby galaxies whereas it is 
not for the Lyman break galaxies). Before trying to quan- 
tify this difference between the escape fractions, we have 
first examined a quantity that is directly observed in both 
cases, the flux ratio /(1500)//(900). 



4.1. The /(1500)//(900) flux ratio of nearby galaxies 

We have calculated the /(1500)//(900) flux ratios in Mrk 
54 and the four nearby star-forming galaxies observed with 
HUT (Table 1). For the comparison with the high-z galax- 
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ies to be meaningful we have to use the /(900) upper limit 
fluxes corrected for gas-phase absorption and foreground 
dust extinction. As for Mrk 54, the latter correction has 
been based on the E{B — V) of Schlegel et al. (199S) given 



in the NED. These corrections are slightly different from 
those of Leitherer et al. (1995) but the conclusions are 
not changed. All steps of calculations are detailed in the 
notes to Table 1. The /(1500)//(900) flux ratios of the 
nearby objects (Table 1 column 7) are found to be larger 
than the value obtained by Steidel et al. (2001) from their 
composite spectrum of Lyman-break galaxies (note that 
their value of 4.6 translates into a value of 1.7 in the units 
of Table 1). This is especially true for Mrk 54 essentially 
because no correction for H I absorption has to be applied. 

Although the LyC photons escaping is probably highly 
variable from galaxy to galaxy, the number of nearby ob- 
jects is sufficient to suggest that the escape is easier at 
redshift « 3.4 than at low redshift. This trend is now sup- 
ported by similarly large /( 1500)// (900) flux ratios re- 



cently reported by Ferguson (2001) for seven galaxies at z 



w 1 in the HDF. The UV luminosity may also be a factor 
in the sense that higher luminosity at w 1500 A would im- 
ply more photons for ionizing the gas and possibly easier 
escape. The nearby galaxies have UV luminosities ranging 
from 7.2 x 10 28 erg s" 1 Hz" 1 at ~ 1500 A for Mrk 54 to 
ten times less. However this largest luminosity is s imilar 
to the demagnified luminosity of cB58 (Pettini et al. 2000 ) 
reported by Steidel et al. as representative of the galaxies 
chosen for the composite spectrum. 



4.2. The /(1500)//(900) flux ratio of Ly break 
galaxies 

We have also examined the /(1500)//(900) flux ratio ob- 
tained for the Lyman break galaxies by Steidel et al. af- 
ter correction for the IGM opacity. This ratio implies an 
escape fraction (when normalized to 1500 A) > 0.5. A 
smaller escape fraction in the absolute sense might be ob- 
tained if only a fraction of the 1500 A photons escape 
without being absorbed by dust. This argument is given by 
Steidel et al. and used by Loeb & Barkana ( 200C ) to infer 
an escape fraction of ~ 10%, close to the limits obtained 
for nearby galaxies. There is however a serious limitation 
with this explanation: a significant dust extinction at 1500 
A would likely imply a larger dust extinction at 900 A and 
in turn a smaller intrinsic /(1500)//(900) ratio, uncom- 
fortably small with respect to those given by stellar pop- 
ulation synthesis models (e.g. Bruzual & Chariot 1993, 
Chariot |l996j , Leitherer et al. |1999| ). The flatness (in f„) 



of the composite spectrum obtained by Steidel et al. also 
suggests a low extinction according to the correlation be- 
tween the ultraviolet extinctio n and the slope [3 of the far- 
UV spectrum (Calzetti et al. |l994[) . The result of Steidel 
et al. is puzzling in another respect. If the escape of LyC 
photons is anisot ropic as discussed above and suggested 
by Pettini et al. (2000), it would be difficult to have only 
the cases with favorable directions appear in a compos- 



ite spectrum. It is also true that the galaxies contributing 
to the composite spectrum have been selected using cri- 
teria that would favor the least extinguished and perhaps 
lowest covering fraction. 

4.3. Quantifying the difference between escape 
fractions 

The LyC escape fraction appears in equation (2) as the 
7(900) flux divided by the Ha flux and is defined by Steidel 
et al. as the /(900) flux normalized by the fraction of 1500 
A photons that escapes. It is tempting to establish a rela- 
tion between these two forms of the escape fraction, inde- 
pendently of the /(900) flux. Equation (2) can be written 
as (neglecting the factors /l and /2) 



fei 



na /(900) /(1500)~ 
°/(1500) /( Ha ) 



On the other hand the escape fraction as defined by Steidel 
et al. is 

/(900) 



csc2 



l.i- 



7(1500) 

with 1.1 being the /(1500)//(900) ratio (in units of Table 
1) as expected from stellar synthesis models and adopted 
by Steidel et al. From the ratio of the two equations we 
get 



csc2 



0.04 



/(Ha) 



2Q0.4A(Hct) 



./;.,• ' /(i5oo)" 

The ratio of the two escape fractions appears as depend- 
ing on the Ha to far-UV (1500 A) ratio, a quantity more 
currently measured than the flux at 900 A. The latter ra- 
tio (expressed in A with an UV flux in erg cm~ 2 s _1 
A -1 ) can take in principle any value from near zero (a 
burst older than 7 Myrs without any ionizing star) to very 
large (dust extinction). However, it is currently measured 
in the range 10 - 40 A (e.g. Buat et al. |1987| , Meure r et 
al. |1999| , Sullivan et al. ^000| , Bell & Kennicutt j200l| ) for 
the bright, low-extinction and UV-selected star- forming 
galaxies as we are dealing with in this paper. The cali- 



brations of Kennicutt (1998) for the far-UV and Ha star 



formation rates correspond to an Ha to far-UV ratio of 
13 A and the model calculations of Glazebrook et al. 
( |1999|) encompass the range 6 - 17 A (extinction not 
included). The / CS c2//csc ratio is therefore found in the 
range 0.4 - 5, assum ing a n extinction at Ha in the range 
0-1.2 (Kennicutt |l998 ). This ratio can be calculated 
for the nearby galaxies and indeed falls into the expected 
range. The large range obtained shows the difficulty of 
comparing escape fractions with different definitions. 



In conclusion, the difference between the large LyC es- 
cape fraction reported at high redshift and the current 
upper limits in nearby galaxies must be interpreted with 
caution because of the difference of the definitions of the 
escape fractions used. In contrast the /(1500)//(900) flux 
ratios in nearby star-forming galaxies suggest that the 
LyC photons escape, at least in the direction of the ob- 
server, less easily from nearby galaxies than from high 
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Table 1. Calculation of the /(1500)//(900) ratio in nearby galaxies (the ratio of 4.6 reported by Steidel et al. is 1.7 
in the units of the table) 



Object 


/(900) 




E{B-V) 


/c(900) 




/(1500) 


/ c (1500) 


/c(1500)// c (900) 




(2) 




(3) 


(4) 




(5) 


(6) 


(7) 


IRAS08339+6517 


< 1.7 x 10" 


15 


0.092 


< 8.8 x 10" 


15 


2.8 x 10" 14 


5.6 x 10" 14 


> 6.4 


Mrk 1267 


< 1.4 x 10" 


15 


0.034 


< 2.6 x 10" 


15 


8.5 x 10" 15 


1.1 x 10" 14 


> 4.2 


Mrk 66 


< 1.6 x 10" 


15 


0.012 


< 2.0 x 10" 


15 


7.5 x 10" 15 


8.2 x 10" 15 


> 4.1 


Mrk 496 


< 1.6 x 10" 


15 


0.020 


< 2.3 x 10" 


15 


8.0 x 10" 15 


9.3 x 10~ 15 


> 4.0 


Mrk 54 


< 4.7 x 10" 


16 


0.015 


< 6.2 x 10" 


16 


2.2 x 10" 14 


2.5 x 10" 14 


> 40 



Notes: all fluxes are in erg cm - s _1 A - . 

Col(2) Corrected for gas-phase absorption (Table 1 , column 12 of Hurwitz et al. 1997 for the galaxies observed with HUT). 
Col(3) Foreground extinction from the NASA Extragalactic Database. Col(4) Corrected for foreground e xtinct ion with 
A — 19.36i?(-B — V), resulting from the extrap olatio n to ~ 940 A of t he extinction law of Cardelli et al. (1989). Col(5) 
From the HUT or IUE spectra (Leitherer et al. 



A = 8.26E{B - V). 



1995 



Kinney et al. 1993). Col(6) Corrected for foreground extinction with 



redshift galaxies. More and better measurements of the 
Lyman continuum of nearby galaxies are needed to dis- 
tinguish a possible evolution effect from selection effects. 

5. Implications 

The stellar contribution to the hydrogen-ionizing back- 
ground has been investigated by a number of authors and 
its intensity anticipated as a function of the LyC escape 
fraction. These results can therefore be briefly revisited in 
the light of present data and especially the /( 1500)// (900) 
ratio. 

W ith calculatio ns sim ilar to those of Madau et al. 
( |1999| ) and Madau ( |2000| ), Steidel et al. directly use their 
measured /(1500)//(900) ratio to translate the luminos- 
ity function of Lyman break galaxies (at 1500 A ) into a 
distribution of Lyman continuum luminosities. They find 
that the 900-A emissivity of galaxies at z ~ 3 exceeds 
the contribution of QSOs by a factor of 5. This result sup- 
ports the role of star formation in the early re-ionisation of 
the hydrogen but relies on the assumption that the mea- 
sured /(1500)//(900) ratio is characteristic of the whole 
population of Lyman-break galaxies. 



At low redshifts Giallongo et al. ( 1997 ) and Shull et al. 
( |l999| ) have shown that the galaxy contribution to the ion- 
izing background rivals that of QSOs for an escape fraction 
of 5%. Such a value is compatible with the present set of 
upper limits obtained with nearby star-forming galaxies. 
In the same vein as above, an upper limit on the 900- 
A emissivity of galaxies at z « can be directly trans- 
lated from the far-UV luminosity density of Sullivan et al. 



( pOOOj ), using the /(150 0)//( 900) ratios in Table 1. As the 
data of Sullivan et al. fl2000| ) are at 2000 A and a mean 
redshift of 0.15, we have first used relations established by 
Cowie et al. ([1999D to get the emissivity 0.9 x 10 26 



erg s 



Hz" 1 Mpc~ 3 at 1500 A and z w (this would be 0.7 x 10 26 
with an (1 + z) 4 evolution factor). With the tightest con- 
straints of Mrk 54 in Table 1 we get a 900-A luminosity 
density < 8 x 10 23 erg s _1 Hz -1 Mpc~ 3 , which is close 
to the emissivity due to galaxies evaluated by Shull et al. 



( |i999l ) for an escape fraction of 5% that would provide a 
comparable background to AGN. 

All the above calculations raise the issue on how the 
value used for the LyC escape fraction is representative 
of the whole population of galaxies. Given all the factors 
(distribution of neutral gas, orientation) that make this 
parameter random, it appears a crude oversimplification 
to use an unique value for all galaxies. At the same time, 
any effort to understand the variations of this parameter 
should carry the same level of difficulty as establishing the 
LyC luminosity function itself. 

As to the present calculations on Mrk 54 and compar- 
ison with other nearby galaxies, it is likely, from the mere 
fact that the LyC escape fraction is random, that any av- 
erage value should lie below the current upper limits. In 
addition, any dust extinction decreases the LyC escape 
fraction below the value expected from neutral hydrogen 
absorption alone (see eq 1 or 2). The nearby galaxies under 
study, selected with significant UV flux for the prospect 
of detecting ionizing radiation, are probably not represen- 
tative of the whole population in terms of dust extinction. 
Such an effect would be present when it comes to using 
an average escape fraction with an Ha luminosity density 
for deriving the LyC luminosity density of galaxies. 



6. Conclusions 

The wavelength domain of the Far Ultraviolet 
Spectroscopic Explorer (FUSE) down to 905 A al- 
lows us to observe the LyC radiation of low redshift 
galaxies above the limit of photoelectric absorption by 
neutral hydrogen in our Galaxy. The high spectral res- 
olution of FUSE combined with an appropriate redshift 
selection allows us to avoid Lyman series absorption by 
neutral hydrogen in our Galaxy over a reasonably large 
wavelength range. The high spectral resolution of FUSE 
is also useful to identify (and eventually to correct for) 
possible gas-phase absorption from molecular hydrogen 
and species other than hydrogen. The observation of the 
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star-forming galaxy Mrk 54 (z = 0.0448) is reported here 
with the following results: 

(i) The LyC flux is not detected. An upper limit can be 
reasonably evaluated by comparison with the background 
in adjacent zones of the detector. It is measured over three 
spectral windows (16 A in total) that are free to a high 
probability of any absorption by neutral hydrogen in our 
Galaxy. No correction for molecular hydrogen absorption 
is needed according to the lack of the corresponding ab- 
sorption features in the spectrum at longer wavelengths. 
A correction for foreground dust extinction is performed. 

(ii) By comparison with the Ha flux of Mrk54, the 
resulting upper limit on the LyC flux translates into an 
upper limit LyC escape fraction / osc < 0.062. This result 
is compatible with current models for the escape of LyC 
radiation in galaxies and adds another significant limit to 
three other cases with limits of the order or less than 10%. 

(iii) A low escape of LyC photons is expected in Mrk 
54 from the large covering fraction of neutral hydrogen re- 
vealed by the Lyman series absorption lines down to Ly7. 
The randomness of the LyC escape fraction as resulting 
from the geometry (anisotropy in escape direction, hole in 
the neutral gas distribution) and large-scale galactic out- 
flows is emphasized. No relation has yet been established 
between the LyC and the Lya emission escape in starburst 
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(iv) The contrast between the detection of aignificant 
LyC flux in a composite spectrum of Lyman break galax- 
ies (Steidel et al. 2001) at z w 3.4 and the fact that only 
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upper limits have been obtained in nearby star-forming 
galaxies is investigated. The /( 1500)// (900) flux ratio is 
found larger in nearby galaxies than in the composite spec- 
trum of Steidel et al. ( |2001 ). The possible trend of LyC 
photons escaping more easily at high redshifts than at 
low redshifts is now supported by the STIS data at z w 
1 of Ferguson ( 2001 ) but more and better observations 
of nearby galaxies are needed to distinguish an evolution 
effect from selection and random effects. 

(v) In addition to an evaluation from the Ha emis- 
sivity and the LyC escape fraction, the 900-A emissivity 
of galaxies and potentially the ionizing background radia- 
tion can now be derived from the measured UV luminosity 
density and (UV to 900 A ) ratio. 

Acknowledgements. We thank the FUSE Team at JHU for 
the suc cessful oporation of euch a complox instrument. J M 



Giallongo, E., Fontana, A., & Madau, P. 1997, MNRAS, 289, 
629 

Glazebrook, K., Blake, C, Economou, F., Lilly, S., & Colless, 

M. 1999, MNRAS, 306, 843 
Gonzalez-Delgado, R. M., Leitherer, C, & Heckman T. 1997, 

ApJ, 489, 601 

Gonzalez-Delgado, R. M., Leitherer, C, Heckman, T. et al. 

1998, ApJ, 495, 698 
Haardt, F., & Madau, P. 1996, ApJ, 461, 20 
Heckman, T. M., Robert, C, Leitherer, C, Garnett, D. R., & 

van der Rydt, F. 1998, ApJ, 503, 646 



Heckman, T. M. 2000, preprint ( |astro-ph/0009075| ) 
Hurwitz, M., Jelinsky, P., & Dixon, W. V. D. 1997, ApJ, 481, 
L31 

Kennicutt, R. C. 1998, ARA&A, 36, 189 

Kinney, A. L., Bohlin, R. C, Calzetti, D., Panagia, N., & Wyse, 

R. F. G. 1993, ApJS, 86, 5 
Kunth, D., Mas-Hesse, J. M., Terlevich, E. et al. 1998, A&A, 

334, 11 

Lockman, F. J., & Savage, B. D. 1995, ApJS, 97, 1 



Loeb, A., fe B arkana, R. 2000, ARA&A pre-print ( |astro- 



D thanks G TCrigg fm- prrwirling gppnifir grrrptg for rpHni-ing 

FUSE data with IRAF, and A. Boselli and G. Gavazzi for 
their optical spectroscopic observation of Mrk 54. This research 
has made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. 



ph/0010467) 



Leitherer, C, Ferguson, H. C, Heckman, T. M., & Lowenthal, 

J. D. 1995, ApJ, 454, L19 
Leitherer, C, Schaerer, D., Goldader, J. D. et al. 1999, ApJS, 

123, 3 

Madau, P., & Shull, J. M. 1996, ApJ, 457, 551 
Madau, P., Haardt, F., & Rees, M. J. 1999, ApJ, 514, 648 



Madau, P. 2000 , Phil. Trans. Roy. Soc. A, in press (astro- 



ph/0003096) 



References 

Barnstedt, J., Gringel, W., Kappelmann, N. 
2000, A&AS, 143, 193 



Meiksin, A., & Madau, P. 1993, ApJ, 412, 34 
Meurer, G. R., Heckman, T. M., & Calzetti, D. 1999, ApJ, 521, 
Grewing, M. 64 

Miralda-Escude, J., & Ostriker, J. P. 1990, ApJ, 350, 1 



Deharveng et al.: Lyman continuum radiation from galaxies 



9 



Moos, H. W., Cash, W. C, Cowie, L. L. et al. 2000, ApJ, 538, 
LI 

Pettini, M., Steidel, C. C, Adelberger, K. L., Dickinson, M., 

& Giavalisco, M. 2000, ApJ, 528, 96 
Savage, B. D., Bohlin, R. C, Drake, J. F., & Budich, W. 1977, 

ApJ, 216, 291 

Sahnow, D. J., Moos, H. W., Ake, T. B. et al. 2000, ApJ, 538, 
L7 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 
5 

Shull, J. M., Roberts, D., Giroux, M. L., Penton, S. V., & 

Fardal, M. A. 1999, AJ, 118, 1450 
Shull, J. M., Tumlinson, J., Jenkins, E. B. et al. 2000, ApJ, 

538, L73 

Songaila, A., Cowie, L. L., & Lilly, S. J., 1990, ApJ, 348, 371 
Stark, A. A., Gammie, C. F., Wilson, R. W. et al. 1992, ApJS, 
79, 77 

Steidel, C. C, Pettini, M., & Adelberger, K. L. 2001, ApJ, 546, 
665 

Sullivan, M., Treyer, M., Ellis, R. S. et al. 2000, MNRAS, 312, 
442 

Tenorio-Tagle, G., Silich, S. A., Kunth, D., Terlevich, E., & 

Terlevich, R. 1999, MNRAS, 309, 332 
Tumlinson, J., Giroux, M. L., Shull, J. M., & Stocke, J. T. 

1999, AJ, 118, 2148 

Vidal-Madjar, A., Kunth, D., Lecavelier des Etangs, A. et al. 

2000, ApJ, 538, L77 



